Asymmetric Wave Propagation in Nonlinear Systems 
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A mechanism for asymmetric (nonreciprocal) wave transmission is presented. As a reference sys- 
tem, we consider a layered nonlinear, non mirror-symmetric model described by the one-dimensional 
Discrete Nonlinear Schrodinger equation with spatially varying coefficients embedded in an other- 
wise linear lattice. We construct a class of exact extended solutions such that waves with the same 
frequency and incident amplitude impinging from left and right directions have very different trans- 
mission coefficients. This effect arises already for the simplest case of two nonlinear layers and 
is associated with the shift of nonlinear resonances. Increasing the number of layers considerably 
increases the complexity of the family of solutions. Finally, numerical simulations of asymmetric 
wavepacket transmission are presented which beautifully display the rectifying effect. 

PACS numbers: 05.45.-a 



The design of small devices to control energy and/or 
mass flows at different spatial scales is a suggestive chal- 
lenge from both a theoretical and applied viewpoint. 
Since the first proposal of a thermal diode pQ , capable of 
transmitting heat asymmetrically between two temper- 
ature sources, several studies appeared in the literature 
[2] , including some experimental realizations [3] [4] . 

A related issue concerns the possibility of devising a 
"wave diode" in which electromagnetic or elastic waves 
are transmitted differently along two opposite propaga- 
tion directions. A so-called optical diode has been pro- 
posed in Ref. [5j |6] and later on discussed both theoreti- 
cally P3 E] and experimentally [9 . There is also a recent 
proposal for a diode based on left-handed metamaterials 
[TO] . Another domain of application is the propagation 
of acoustic pulses in granular systems. Indeed, experi- 
mental studies demonstrated a change of solitary wave 
reflectivity from the interface of two granular media [TT] . 

The basic idea for a "wave diode" was proposed for 
nonlinear photonic crystals [7 : the rectification depends 
on whether the second harmonic of the fundamental wave 
is transmitted or not. As a definition of a rectifying 
device we would instead propose that the transmitted 
power at fixed incident amplitude and at the same fre- 
quency uo should be sensibly different in the two opposite 
propagation directions. In the present Letter we pur- 
sue a novel approach that exploits distinctive features of 
nonlinear dynamical systems such as multistabiliy and 
amplitude-dependent resonances to achieve such an ef- 
fect. 

In a linear, time-reversal symmetric system this possi- 
bility is forbidden by the reciprocity theorem [12] . There- 
fore, one needs to consider nonlinear and asymmetric 
systems [13] . As a reference model we will focus on the 
Discrete Nonlinear Schrodinger (DNLS) equation p~4j[15] 
with spatially varying coefficients. It has been demon- 
strated [16] that DNLS can be a sensible approximation 



for the evolution of longitudinal Bloch waves in layered 
photonic or phononic crystals (Fig. [I]). Variable coeffi- 
cients describe different nonlinear properties of each layer 
and the presence of defects. Beyond its relevance in many 
different physical contexts, the DNLS equation has the 
big advantage of being among the simplest dynamical 
systems amenable to a complete theoretical analysis. For 
our purpose, it is particularly convenient as it allows to 
solve the scattering problem exactly without the compli- 
cations of having to deal with wave harmonics. 
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FIG. 1: (Color online) Sketch of a layered photonic or 
phononic system. The central N layers are nonlinear, non 
mirror-symmetric with respect to the structure center. In the 
limit of a vanishing width of the thinner layers the dynamics 
within high-frequency Bloch bands is described by a DNLS 
equation for the envelope ip n (see Refs. [16, 17 for details). 



More precisely, let us consider the stationary DNLS 
equation defined on an infinite one-dimensional lattice 

W&n = V n 1p n ~ - + a n \l/j n \ 2 l/j n . (1) 

We will assume the usual scattering setup where V n and 
a n are nonvanishing only for 1 < n < N . The two semi- 
infinite portions (n < 1, n > N) of the lattice, model two 
leads where the wave can propagate freely. 
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Let us look for solutions of the associated transmission 
problem 
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where uj = — 2cos/c and < k < tt for the wave com- 
ing from the left direction; Rq, R and T are the incident, 
reflected and transmitted amplitudes respectively. The 
solution sought must be complex in order to carry a non- 
vanishing current J = 2\T\ 2 sin/c. 

To break the mirror symmetry with respect to the cen- 
ter of the nonlinear portion, one must choose at least 
one of the two sets of coefficients V ni a n such that 
V n Vn-u+i, a n 7^ OLN-n+i- Note that the transmis- 
sion of the right-incoming wave with the same Ro and 
uj is computed by solving the problem with (V n: a n ) — > 
(V/v_ n+ i, c^jv-n+i) (i- e - "flipping the sample"). In the 
following, we will refer to such solutions as those having 
negative wavenumbers, —k. Nonlinearity is essential as 
for a n = the transmission coefficient is the same for 
waves coming from the left or right side, independently 
on V n due to time-reversal invariance of the underlying 
equations of motion [T8] . 

The standard way to solve the problem is to introduce 
the (backward) transfer map [T71 [T9H2T] 



-v n + 04 - uj + a n \u n \ 2 )u n , 



(3) 



where u n = ifj n and v n = ^n+i- Note that these 
are complex quantities therefore the map is nominally 
four-dimensional. However, due to conservation of en- 
ergy and norm, it can be reduced to a two-dimensional 
area-preserving map [T7J fT9H2T] with an additional con- 
trol parameter (the conserved current J). The solu- 
tions are found by iterating (|3| from the initial point 
u N = Tex.p(ikN), v N = T exp(ik(N + 1)) dictated by 
the boundary conditions of Eq. For fixed T and /c, 
the incident and reflected amplitudes are determined as 



Ro = 



ex.p(—ik)uo — vo 
exp(— ik) — exp(z/c) ' 



R = 



ex.-p(ik)uo — vo 
exp(z/c) — exp(— ik) 



and the transmission coefficient is £(&, |T| 2 ) = |T| 2 /|i?o| 2 
Note that if (iio,^o) = (un^n) (periodic point of the 
map) then t = 1. 

In order to illustrate the effect in the simplest case we 
consider the DNLS dimer N = 2. The coefficient t can be 
thereby computed analytically iterating the map twice: 



-ik 



1 + (y - e ik )(e ik - 5) 



(4) 



where 6 = V 2 -uJ + a 2 T 2 , v = V x -cj + aiT 2 [l-2£cosfc + 
S 2 ]. The formulas applies for k > 0. As explained above, 
to solve for the case k < one has to exchange the sub- 
scripts 1 and 2 which is equivalent to reverse the sample. 

Up to now no hypothesis has been made on the cou- 
pling coefficients. For simplicity, we impose = ct > 



henceforth and let Vi } 2 = Vb(l ± e). I n the symmetric 
case e = 0, it can be verified that transmission is unity 
for V + otT 2 = (for V < 0) and V + aT 2 = uj (for 
Vo < uj). These nonlinear resonances can be regarded 
as the continuation of the extended (i.e. perfectly trans- 
mitting) states of the corresponding linear problem [22] 
to nonvanishing a. As a result, the transmission curves 
(I^Rol 2 , 1^1 2 ) are tangent to the bisectrix in two points 
where t = 1 (dashed lines in Fig. ^ and b respectively). 

For the asymmetric case, e ^ 0, the resonances are 
detuned differently for the k > and k < cases, lead- 
ing to a nonreciprocal transmission (see the solid lines 
in Fig. [2^i,b) with maximal t smaller than one. In par- 
ticular, there exist two windows Wi^ of amplitude of 
the incident wave Ro where we have three solutions for 
k > and only one for k < 0. Here the asymmetry 
become maximal. Fig. [3] illustrates the different types 
of solutions in the region W\ (see the points in the in- 
set of Figj2^i). The two low-amplitude solutions are very 
similar as expected for weak nonlinearity [23] . 
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FIG. 2: (Color online) DNLS dimer, N = 2, V = -2.5, a = 1, 
\k\ = 0.1. Comparison between the symmetric case (s = 0, 
dashed lines) and the asymmetric one (e = 0.05, solid lines), 
(a) Transmission curves, the inset is an enlargement of the 
low-amplitude region. The vertical lines mark the turning 
points of the curves. Accordingly, the heavy lines on the 
horizontal axis are the multistability windows Wi,2 where the 
diode effect maximally occurs, (b) Transmission coefficients 
as a function of the transmitted intensity. 



Since the phenomenon is of nonlinear origin the asym- 
metry depends on both frequency and amplitude. To 
quantify its efficiency, in Fig. [4] we report the rectifying 
factor 



/ = 



t (k, in 2 )-t(-fc,iri 2 ) 
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(5) 
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FIG. 3: (Color online) Square modulus of the solutions corre- 
sponding to the same incident amplitude \Ro\ 2 = 0.4 (marked 
by full dots in the inset of Fig.[2|, other parameters as in previ- 
ous figure. Upper panel: the three left-propagating solutions 
corresponding to |T| 2 = 0.327,0.377,0.01 (top to bottom). 
Lower panel: the right-propagating solution corresponding to 



0.01. The oscillations are caused by the interference 



between the incident and reflected waves (wavenumber 2k). 



which approaches ±1 for maximal asymmetry. Note that, 
although increasing e broadens the regions in which |/| 
is relatively large, the overall transmitted intensity is re- 
duced as well. 
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FIG. 4: (Color online) Contour plot of the rectifying factor 
|5| for increasing asymmetry level e. Other parameters as in 
the previous figure. 



Increasing the number N of nonlinear layers will con- 
siderably increase the complexity of the transmission pat- 
terns. Indeed, for large TV, the transmitting (resp. non- 
transmitting) solutions correspond to bounded (resp. es- 
caping) orbits of the map |3|. Due to the mixed phase- 
space, t is nonvanishing only on a fractal set of the param- 
eter space [19]. The system thus exhibits a complicated 
form of multistability, meaning that for a fixed input one 



may have (infinitely) many outputs jTTl [19] . We thus ex- 
pect that the structure of the windows W n will approach 
that of a Cantor set in the large N limit. Indeed, Fig. [5] 
indicates that even for small iV, the W n shrink and the 
texture of the regions where |/| < 1 complicates consid- 
erably. 
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FIG. 5: (Color online) Upper panel: transmission curves for 
DNLS with N = 4, linearly-modulated parameter V n — Vb[l+ 
e(l — 2(n — 1)/(N — 1)] , other parameters as in the previous 
figures. Upper panel: comparison between the symmetric case 
(e = 0, dashed lines) and the asymmetric one (e — 0.05, solid 
lines). Lower panels: contour plots of the rectifying factor 
([5| for increasing asymmetry. Similar results are obtained for 
different types of modulation. 



What are the consequences of the above results on 
the transmission of wavepackets? In a nonlinear system 
where the superposition principle no longer holds, the 
connection between the two problems is not trivial. To 
address this problem, we solved numerically the time- 
dependent DNLS 

l(j) n V n <t> n ~ 0n+l - 0n-l + OLn\<t>n\ 2 (t>n (6) 

on a finite lattice \n\ < M with open boundary condi- 
tions, for the case of the dimer discussed above. We take 
as initial condition a Gaussian wavepacket 

n ) 2 
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(7) 



The upper panels of Fig. [6] display the evolution of two 
packets with the same I and opposite wavenumber ko 
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impinging on the nonlinear dimer. The asymmetry of 
their propagation is manifest. Here, the parameters have 
been chosen empirically to obtain the maximal asymme- 
try. More precisely, we first measured the wavepacket 
transmission coefficient t p , defined as the ratio between 
the transmitted norm J2 n >N at the enc ^ °f the run 
divided by the initial one X^ n <o 1^(0) | 2 , as a function 
of \I\ 2 . The data of Fig. [6] correspond to the ampli- 
tude / for which the transmission is maximal for the 
left-incoming packet (t p ~ 0.8) and minimal for the right- 
incoming one (t p c± 0.3). Although the packets are sig- 
nificantly distorted after scattering, the Fourier analysis 
shows that they remain almost monochromatic at the in- 
cident wavenumber ko (lower panels of Fig. [6| . It is also 
noteworthy that reflection is associated with the creation 
of a localized excitation, strongly reminiscent of the non- 
linear impurity modes [T5] . 

In conclusion, based on the DNLS model, we have 
demonstrated a mechanism which leads to nonreciprocal 
wave transmission. The new class of solutions found here 
are of interest both theoretically as well as to envisage 
possible experimental realizations in nonlinear layered 
photonic or phononic systems. To the extent to which 
the DNLS can be considered a realistic model for such 
media, our results may open the way to novel strategies 
to control and optimize wave propagation and to design 
devices for sound or light rectification. 
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FIG. 6: (Color online) Numerical simulations of the propaga- 
tion of Gaussian wavepackets, Eq. Q impinging on a DNLS 
dimer. Here V -2.5, \k \ 1.57, s 0.05, M 500, 
\I\ 2 — 3, w — 10 4 and no = T250, respectively. Lower panels: 
Power spectra of the real part of n at times t — (green 
dashed line) and t = 250 (red solid line). 
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